Vibrational wave packet induced oscillations in two-dimensional electronic spectra. 

I. Experiments 
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I. INTRODUCTION 

Coherent wave packet generation in the manifolds of 
molecular exited states takes place if the excitation pulses 
are spectrally broad enough to coherently excite a num- 
ber of vibrational or electronic levels. Femtosecond (fs) 
laser pulses meet this criterion in a number of molecu- 
lar systems and experiments employing such pulses often 
show signatures of wave packet motion. Modulations of 
signal intensities induced by vibrational wave packet mo- 
tion have been observed in a number of nonlinear spec- 
troscopic experiments, including pump- probe (PP) [l-l3|, 
transient grating (TG) [4] , and a variety of photon-echo 
spectroscopies [J-lsl- A number of theoretical investiga- 
tions supported the conclusions drawn in these exper- 
iments |9l-[l3fl. The influence of intramolecular vibra- 
tions should be expected consequently also in the recently 
implemented heterodyne-detected two-dimensional elec- 
tronic spectroscopy (2D-ES) technique (l3 - [2H . which is 
based on the correlation of electronic coherences evolving 
during two time periods. 2D electronic spectra recorded 
for distinct waiting times provide information on coupling 
patterns and evolution of populations and coherences via 
line-shape modulations and the appearance and disap- 
pearance of cross peaks [H [11, l20l i22l425t . 

Recent theoretical works on model dimer [2^ and 
trimer [28j systems, and experimental investigations on 
excitonic complexes [l^HOl and conjugated polymers (3l| 
reported on signatures of electronic wave packets in 2D 



electronic spectra resulting from coherent excitation of 
several electronic levels. Possible enhancement of the en- 
ergy transfer efficiency in photosynthetic light harvesting 
complexes, due to the presence of excitonic wave pack- 
ets, was suggested by Engel et al. [2^ . An assignment of 
the observed spectral modulations to an excitonic wave 
packet motion was made, based on a good agreement of 
the experiment with the predictions of an excitonic model 
[26l.[29|. In the general case, however, contributions of vi- 
brational wave packet motion cannot be excluded and it 
is therefore of high importance to study the vibrational 
case separately. To the best of our knowledge, vibra- 
tional effects in 2D-ES have been addressed so far only 
theoretically [s^, 'ss'l and in our recent work on the same 
molecular system [20]. 

To study the effect of vibrational coherence with- 
out perturbations from additional electronic levels, a 
two-level system is preferential. Quantum-chemistry 
calculations (cf. Section II) indicate that N,N'-bis(2,6- 
dimethylphenyl)perylene-3,4,9,10-tetracarboxylicdiimide 
(PERY), when excited at 20000 cm~*, can be approx- 
imated as a two-level electronic system with a number 
of vibrational modes coupled to the electronic transition 
[3^ . In the linear absorption spectrum of PERY in 
solution four well resolved vibrational peaks with a 
frequency separation of «1400 cm~* can be discerned 
(Fig. Uti-l- contrast to a recent work by Tekavec 
et al. [35|, in which white light pulses where used 
to probe the vibronic dynamics, this progression can 
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Figure 1: (a) Linear absorption spectrum of PERY in toluene 
(black solid line) and spectrum of the excitation pulses (green 
dashed line). The red line indicates the position of the low- 
est energy transition obtained from ZINDO/S calculations. 
The inset shows the chemical structure of PERY. (b) Huang- 
Rhys factors determined from semiemprical AMI calcula- 
tions. Only modes that posses Huang-Rhys factors larger 
than 0.005 are plotted. 



not be covered with the bandwidth of our laser pulses 
(900 cm~^). However, by tuning the central excitation 
energy approximately to the center of the lowest energy 
peak, nonlinear spectroscopic signals such as transient 
grating (TG) and three-pulse photon echo peak shift 
(3PEPS) have been shown to be modulated by a low 
frequency (140 cm~^) vibrational mode (s^ . 

In this contribution we investigate experimentally how 
coherent excitation of several vibrational levels influences 
2D electronic correlation and relaxation spectra. To this 
end 2D electronic spectra of PERY have been recorded 
for 20 waiting times t2 ranging from to 800 fs covering 
several periods of pronounced vibrational modulations of 
the spectra. In a recent work on the same molecular sys- 
tem we have shown that this low frequency mode induces 
a periodic change in the ellipticity of the real (absorptive) 
part and in the slope of the nodal line separating the 
positive and negative contribution of the imaginary (dis- 
persive) part [20]. Building up on these qualitative ob- 
servations, in this contribution we will assign these mod- 
ulations to periodic enhancement and decrement of the 
rephasing and non-rephasing signal parts and establish 
connections to the frequency-frequency correlation func- 
tion M(t). In the accompanying paper Ref. [36] (paper 



II) the underlying theory is presented. 

This paper is organized as follows. We start in Sec- 
tion II with a quantum- chemical characterization of the 
molecule under investigation. The experimental details 
are provided in Section III. In Section IV two-dimensional 
electronic spectra are presented and connections to the 
frequency-frequency correlation function are established. 
We close in Section V with some concluding remarks. 
The appendix contains a comparison to related one- 
dimensional four- wave mixing signals. 



II. QUANTUM-CHEMICAL CALCULATIONS 

This section characterizes PERY theoretically and 
spectroscopically with respect to its two-dimensional tt- 
conjugation by means of quantum-chemical calculations. 
The electronic ground state geometry of PERY has been 
optimized using Density Functional Theory (DFT) fs^ 
based on the Becke's three parameter hybrid functional 
using the Lee, Yang and Parr correlation functional 
(B3LYP) [H]. The polarized split-valence SV(P) basis 
set [3^ has been used. The obtained optimal structures 
were checked by normal mode analysis (no imaginary 
frequencies for all optimal geometries). Based on the 
optimized geometry, the vertical transition energies and 
oscillator strengths between the initial and final states 
have been calculated us ing the Configuration Interaction 
Singles (CIS) method [40| based on the semiempirical 
ZINDO/S (Zerner's Intermediate Neglect of Differential 
Overlap) Hamihonian [U. For the ZINDO/S calcula- 
tions, the single excitations from the 10 highest occupied 
to the 10 lowest unoccupied molecular orbitals were con- 
sidered. The DFT calculations were done using the Tur- 
bomole 5.7 package [42] and the optical transitions were 
computed using the Argus Lab software [i^. The evalu- 
ated theoretical values are used for the interpretation of 
experimental spectroscopic measurements. 

The optimal geometry of D2 symmetry exhibits a pla- 
nar central part (Fig. [2^). The orientation of the lat- 
eral aromatic rings and the core skeleton is nearly per- 
pendicular (dihedral angles of 81.2°) causing a negligible 
TT-conjugation between these molecular fragments. The 
lateral rings are sterically hindered due to the presence 
of methyl groups in ortho-position. The opposite mutual 
arrangement of the neighboring lateral rings is not par- 
allel as can be perceived from the right part of Fig. 

The photophysical properties of PERY are primarily 
determined by the planar geometry of the central part 
and the related size of the 7r-conjugation. The experi- 
mental absorption spectrum of PERY in toluene (Fig.[T^) 
can be characterized by an intense absorption band with 
four peaks which are detected in the spectral range of 
18000 to 24000 cm^i. The equidistant energy differ- 
ence (1400 cm~^) between the four peaks indicates a 
harmonic character of the vibrational progression. The 
calculated vertical ZINDO/S excitation energies and the 
oscillator strengths are in good agreement with the ex- 
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Figure 2: (a) Side and front views on the B3LYP/SV(P) op- 
timal structure of PERY in the electronic ground state, (b) 
Visualization of the HOMO (bottom) and LUMO (top) ob- 
tained from ZINDO/S//B3LYP/SV(P) calculations. 



perimental absorption spectrum. The ZINDO /S method 
was parametrized for spectroscopic properties of the tt- 
conjugated aromatic systems and gives transition ener- 
gies of 18725 cm~^ with an oscillator strength f= 1.7814. 
28631 cm-i (f= 0.0001), 31929 cm~^ (f= 0.0002), and 
34536 cni^^ (f= 0.8228). The lowest energy transition is 
depicted as red line in Fig. [1^. The absence of other in- 
tensive electronic transitions in the region between 18000 
and 24000 cm~^ indicates the vibronic origin of the pro- 
gression observed in the experimental absorption spec- 
trum. Therefore, the first and most intensive peak in the 
experimental record can be attributed to the vi- 
bronic transition. In this context it is useful to examine 
the ZINDO/S HOMO and LUMO of the molecule. The 
HOMO-LUMO excitation plays a dominant role in this 
energetically lowest optical transition according to the 
frontier molecular orbital theory, with its contribution 
being more than 30%. These selected ZINDO/S orbitals 
are visualized in Fig. demonstrating that both or- 
bitals are delocalized over the center of the core part only. 
The lobes of the LUMO are oriented more or less per- 
pendicular to the lobes of the HOMO. The shape of the 
LUMO compared to those of the HOMO indicates that 
the optical excitation is spread from the central atoms of 
the core towards the nitrogen atoms. 

In order to determine the origin of the vibronic motion 
reflected in the linear and nonlinear spectroscopic experi- 
ments, vibrational line spectra have been calculated. The 
complete set of computed vibrational frequencies and in- 
tensities together with a visualisation of the most domi- 



nant modes and Huang-Rhys factors are collected in the 
Supporting Information. For determining the Huang- 
Rhys factors semiemprical AMI calculations were per- 
formed using the Mopac program package 2002 [44*1 . The 
calculated vibrational spectrum possesses 210 fundamen- 
tals between 13 and 3243 cm^^. With the limited band- 
width of our excitation pulses we can only cover modes up 
to frequencies of approximately 900 cm~^. From Fig.[T|3, 
in which modes with a Huang-Rhys factor larger than 
0.005 are plotted, it becomes obvious that only a few 
modes contribute significantly to the observed spectro- 
scopic response. 



III. EXPERIMENTAL 

PERY was purchased from Sigma-Aldrich and used 
without further purification. Solutions with a concen- 
tration of 3 • 10"** M were prepared by dissolving PERY 
in toluene (spectrophotometric grade, Merck), sonication 
for wlO minutes, and filtration to remove undissolved 
residues. 

Pulse generation is accomplished by a regenerative 
titanium-sapphire amplifier system (Reg A 9050, Coher- 
ent Inc.) operating at a repetition rate of 200 kHz. Con- 
version into the visible spectral region is attained with 
a noncollinear optical parametric amplifier (NOPA) [i^ . 
For our present study the central excitation energy was 
set to 18800 cm-i (FWHM 920 cm^^) and the pulses 
were attenuated to 10 nJ per pulse at the sample po- 
sition. Second-order dispersion is eliminated with a se- 
quence of fused-silica prisms, while chirped mirrors com- 
pensate for third-order dispersion, yielding a pulse du- 
ration of 16 fs. The pulse duration is determined with 
an improved version of spectral phase interferometry for 
direct electric field reconstruction (SPIDER), zero addi- 
tional phase (ZAP) SPIDER [Hi, and cross-checked with 
SHG-FROG [47|. 

Our approach for recording 2D electronic spectra is 
based on a diffractive-optics based setup following the 
experimental configuration reported in Ref. [l^ [23, |48| 
(Fig. EJi). A 50/50 beam splitter (BS) splits the NOPA 
output into two beams of equal intensity, one of which 
can be delayed with respect to the other by a computer 
controlled delay stage to introduce delay t2 (cf. Fig. [Hb 
for a definition of time delays). A diffractive optical el- 
ement (DOE), optimized for diffraction into ibfirst or- 
der, is used to split the two beams focused onto it into 
four beams fl^. From this point on, all four beams are 
reflected from common mirrors to maintain the phase 
correlations jH, [l3|- A spherical mirror (SM) after the 
transmission grating parallelizes the four beams which 
are arranged to form the four corners of a square. Pulse 
k2 passes through a pair of glass wedges (WP) that are 
oriented in an anti-parallel fashion |17| . By moving one 
of the wedges with a computer controlled delay stage, 
the amount of glass in the beam path and therefore the 
delay of the pulse is changed without introducing lateral 
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Figure 3: (a) Schematic setup for recording 2D electronic 
spectra. BS . . . beamsplitter, CP . . . compensation plate, SM 
. . . spherical mirror, DOE . . . difTractive optical element, WP 
. . . wedge pair, ND . . . neutral density filter, S . . . sample, L . . . 
lens, CCD . . . charge coupled device camera, (b) Designation 
of pulses and time delays in four wave mixing experiments. 



beam shifts. This configuration allows for precise delay 
movements (5.3 as) without infiuencing the phase stabil- 
ity of the phase locked pulse pair and is used to introduce 
delay ti. Pulses ki and ks pass through equivalent but 
static wedge pairs to balance the dispersion. The LO is 
attenuated by a neutral density filter (ND) by four to 
five orders of magnitude and reaches the sample approxi- 
mately 550 fs before the other pulses. All four beams are 
focused by a spherical mirror onto the sample to generate 
the nonlinear signal. In order to avoid unwanted signals 
from a cuvette window and to prevent precipitation of the 
sample on the cell windows, a wire-guided gravity driven 
jet [49l [soj is implemented, yielding a film of approxi- 
mately 200 /im thickness. Due to the square geometry of 
the excitation beams and the LO, the signal propagation 
direction is coincident with the propagation direction of 
the LO allowing for heterodyne detection of the signal 
field with the LO. Data acquisition is accomplished by 
means of a thermoelectrical cooled charge-coupled-device 
(CCD) spectrometer. A 2400 lines/mm grating in com- 
bination with a 1024 pixel array allows for a spectral 
resolution of 3.3 cm^^ in uj^. 

For recording 2D electronic spectra, t2 is fixed at a 
certain value and ti is scanned over ±100 fs with a step 
size of 0.65 fs to fulfill the Nyquist sampling criterion. 
The electric field of the signal as a function of frequency 
(W3) is reconstructed by Fourier transform of each inter- 
ferogram (for every ti-step), filtering and back Fourier 
transform [5l|. Fourier transform along ti yields the sec- 
ond frequency dimension wi . The frequency resolution in 
oji is solely determined by the maximal temporal separa- 
tion between pulses ki and k2 and amounts to 26 cni~^ 



in our experiments. In a final step the absolute phase of 
the 2D spectrum has to be determined. Since the PP sig- 
nal does not depend on constant phase shifts of either the 
pump or the probe beam, it can be used for this purpose 
[52]. To this end the projection slice theorem is applied, 
that states that the projection of the real part of the 
2D spectrum onto ^3 equals the spectrally resolved (SR) 
PP signal. By multiplying this projection by a factor of 
e"*"^ and adjusting cj) so that the projection matches the 
SRPP signal, the absolute phase of the 2D spectrum is 
recovered. 



IV. RESULTS AND DISCUSSION 
A. Two-dimensional electronic spectra 

With second-order nonlinear signals being symmetry- 
forbidden in isotropic media such as bulk liquids, the 
third-order nonlinearity is the lowest nonlinear order at 
which system-bath interactions and solute dynamics in 
liquids can be probed. The third-order nonlinear polar- 
ization pf^^^ (ii, ^2, ^3) is thereby induced by three inter- 
actions with pulsed electromagnetic fields (with wavevec- 
tors ki, k2, ks and temporal separation ti and ^2) and 
it radiates the signal in the phase-matched directions 
±ki ± k2 ± ks. Depending on the scanning procedure, 
the detection direction (phase matching), and the detec- 
tion scheme (time vs. frequency resolved; homodyne vs. 
heterodyne detection) , different techniques can be distin- 
guished. In this subsection we focus on two-dimensional 
electronic spectroscopy, related one-dimensional four- 
wave mixing techniques are discussed in the Appendix. 

In contrast to one-dimensional FWM techniques, 
where time- or frequency-integrating detection schemes 
project the full information content onto a single axis (i.e. 
intensity vs. time or frequency), 2D spectra provide a 
more complete picture by extending the dimensionality of 
representation. Selective excitation, competing with high 
temporal resolution in one-dimensional experiments, is 
not an issue in 2D spectroscopy owing to the mathemat- 
ical properties of Fourier transform, i.e. the frequency 
resolution being determined by the scanned time range 
and vice versa. 2D electronic spectra are recorded by 
scanning ti over a defined interval for a fixed value of ^2 
and recording the signal in the phase matched directions 
±ki =p k2 + ks in a heterodyne detection scheme. Het- 
erodyne detection serves not only to boost weak signals, 
but also to extract the amplitude and phase of the signal 
field by Fourier transform algorithms [l^ . [5l|. Fourier 
transform with respect to ti and ^3 results in a complex 
signal of the form 

Jo Jo 

X e+''^^*3e±^cuiti^^^^^^^ (^^^ 
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where the negative sign in the exponent accounts for the 
rephasing signal part (i.e. ti > 0, kg = — ki + k2 + ks) 
and the positive sign for the non-rephasing signal part 
(i.e. ti < 0, kg = +ki — k2 + ka). Complex signals can 
be represented either by their real and imaginary parts or 
by their amplitude and phase. These two representation 
schemes are displayed for PERY in Figs. S] and [S] for t2- 
times of 100, 200, 300, 450, 550, 650, and 800 fs. 

The real part of a 2D spectrum (first column in Fig.|4^) 
is the product of absorptive line-shapes in lui and in 
UJ3. In agreement with approximating PERY as a two- 
level electronic system, where ground state bleaching and 
stimulated emission are the only two contributions, it 
shows a single positive feature. The maximum of the 
peak is centered slightly below the diagonal (|wi| > wa), 
indicative of a small Stokes shift. The peak experi- 
ences periodic modulations in its shape with a period 
of 240 fs (a detailed analysis is presented in Fig. [7]). For 
i2-times of 200, 450, and 650 fs, the shape of the ab- 
sorptive peak is strongly elliptical with the major axis 
of the ellipse oriented along the diagonal. Contrarily, 
for t2 = 100, 300, 550, and 800 fs the absorptive peak 
adopts a more circular shape and the major axis is ori- 
ented nearly parallel to the oja-axis. 

The imaginary part of a 2D spectrum is the product of 
an absorptive line-shape in wi and a dispersive line-shape 
in 0^3. In case of a two- level system the imaginary part 
features two contributions of opposite sign (first column 
in Fig. [3)3). Similar to the real part, the imaginary part 
too exhibits oscillations with a period of 240 fs. The most 
obvious oscillation is seen in the nodal line separating the 
positive and negative feature. The angle of this nodal line 
with respect to the wi-axis changes from being positive 
for t2 — 200, 450, and 650 fs, to zero or even negative 
values for ^2 = 100, 300, 550, and 800 fs (cf. Fig. [7] for a 
detailed evaluation) . In addition to this oscillation of the 
nodal line, the width of both contributions periodically 
changes from being narrow for ^2 = 200, 450, and 650 fs, 
to being broad for ^2 = 100, 300, 550, and 800 fs. 

In Fig. [S] we show amplitude and phase 2D spectra of 
PERY for i2-delays of 100, 200, 300, 450, 550, 650, and 
800 fs. The amplitude spectra (first column in Fig. [S^) 
resemble the real part spectra shown in Fig. 2^, with a 
single feature centered slightly below the diagonal. Sim- 
ilar to the real part spectra this positive peak features 
oscillations in its orientation and shape. For ^2 = 200, 
450, and 650 fs the peak shape is highly elliptical and its 
major axis is oriented along the diagonal. For t2 = 100, 
300, 550, and 800 fs the peak acquires a more circular 
shape and is oriented nearly parallel to the wa-axis. 

In the phase spectra (first column in Fig. [SJj) the lines 
of constant phase exhibit oscillations in their orientation. 
For t2 — 200, 450, and 650 fs the lines of constant phase 
are oriented along the diagonal, whereas they turn to be- 
come horizontal at t2 — 100, 300, 550, and 800 fs. Note 
that the phase spectra are shown for a narrower range 
than the real, imaginary, and amplitude part spectra. All 
the observed line-shape modulations will be assigned be- 



low to periodic modulations in the strength of the rephas- 
ing and non-rephasing contributions to the total signal. 



B. Rephasing and non-rephasing signal parts 

The third-order nonlinear polarization -P*-'^^(ti, ^2, ^3), 
which underlies all FWM signals, can be expressed 
as a convolution of the appropriate response functions 
^i(^i,i2ji3) with the electric fields of the excitation 
pulses |53| . The four Feynman diagrams graphically il- 
lustrating the four response functions for a two-level sys- 
tem within the rotating wave approximation are depicted 
in Fig. [6l In these diagrams the left and right vertical 
lines indicate the evolution of the \ket) and {bra\ of the 
density matrix, respectively. Wavy arrows denote inter- 
actions with the electromagnetic field, and time evolves 
from bottom to top. In all four diagrams, the first in- 
teraction of the laser field with the molecular system in- 
duces a coherence between the ground state (g) and the 
first excited state (e) in which the system evolves for a 
time ti. Interaction with the second laser pulse creates a 
population or vibrational coherence in either the ground 
(diagrams R3 and R4) or the excited (diagrams Ri and 
R2) state. After a time t2 the third pulse again induces 
a coherence between the ground and excited state that 
radiates the signal field in the phase-matched direction 
kg. 

The four Feynman-diagrams shown in Fig. [5] can be 
classified as ground state bleaching (GSB) and stimu- 
lated emission (SE), depending on whether the system 
evolves in the ground state (i?3, R4) or in the excited 
state (i?i, R2) during t2, respectively. A different clas- 
sification of these four diagrams arises from the ordering 
of the first two interactions. In diagrams R2 and Rt, the 
first interaction takes place with pulse — ki, whereas in 
diagrams Ri and R4 the first interaction takes place with 
positive sign (-|-ki). As a consequence the system evolves 
in conjugate frequencies during ti (uJge) and (uJeg) in 
diagrams R2 and R3. This allows for rephasing in an 
inhomogeneously broadened ensemble, resulting in the 
formation of an "echo" in case the system has kept some 
memory of its transition frequency over t2. These dia- 
grams therefore contribute to the rephasing signal part. 
In contrast, the system evolves with the same frequency 
during ti and (tdeg) in diagrams i?i and R^. These two 
diagrams are denoted non-rephasing and the systems's 
response is close to a free induction decay. Rephasing 
diagrams contribute to the signal for positive delays ti 
(i.e. pulse 1 preceding pulse 2), whereas non-rephasing 
diagrams contribute for negative delays of ii (i.e. pulse 
2 preceding pulse 1). In addition to these four diagrams 
there are diagrams that contribute to the signal only if 
the pulse separation is shorter than the pulse duration, 
i.e. when the excitation pulses overlap. These diagrams 
will not be taken into account in our discussion. 

Scanning ti over positive and negative values equally 
weights rephasing and non-rephasing pathways. This 
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Figure 4: (a) Real and (b) imaginary part of the 2D spectra of PERY in toluene for t2-delays of 100, 200, 300, 450, 550, 650, 
and 800 fs. The first, second, and third column of each panel display the corresponding total signal, the rephasing, and the 
non-rephasing part, respectively. All total spectra are normalized to their absolute maximum value, whereas the rephasing and 
non-rephasing parts are plotted with their respective contribution to the total signal. Contour lines are drawn at 5% intervals 
starting at ±10%. Red lines indicate positive signals, blue lines negative ones. The solid line in the real part indicates the 
diagonal, whereas the black lines in the imaginary part are drawn at the zero crossings between the positive and negative 
features. 



scanning procedure has been shown to produce purely ab- 
sorptive real parts and purely dispersive imaginary parts, 
that are linked via the Kramers-Kronig relation jsj, ■ 
The rephasing and non-rephasing spectra on their own 
exhibit phase-twisted line-shapes resulting from mixing 



of dispersive and absorptive features. The dissection of 
the total 2D spectrum into its corresponding rephasing 
and non-rephasing parts is shown in the second and third 
row of Figs. [Hand [5] for the real part, imaginary part, am- 
plitude, and phase, respectively. The decomposition was 
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Figure 5: (a) Amplitude and (b) piiase of the 2D spectra of PERY in toluene recorded at f2-delays of 100, 200, 300, 450, 
550, 650, and 800 fs. The first, second, and third column of each panel display the corresponding total signal, the rephasing 
and the non-rephasing part, respectively. The amplitude part total spectra are normalized to their absolute maximum value, 
whereas the rephasing and non-rephasing parts are plotted with their respective contribution to the total signal. Contour lines 
are drawn at 5% intervals starting at ±10% in the amplitude part spectra and at intervals in the phase spectra. The 
zero-phase line is indicated in each phase spectrum. Note the different spectral range in the representation of the phase spectra. 



performed by recording the full scan and evaluating only 
positive ti-values for the rephasing part and only neg- 
ative ii-values for the non-rephasing part (56| . Strictly 
speaking, rephasing signals show up at {— wi, -l-ws} ac- 
cording to the conjugate frequency evolution during t\ 
and t'i. Non-rephasing signals on the other hand ap- 



pear at {-l-wi, -l-wa}, since the system evolves with the 
same frequency during t\ and ^3. However, to facilitate 
comparison of the rephasing, non-rephasing, and total 
spectra, we plot all signals along {— cl>i,+W3} [stI. [58|. 

From the inspection of Fig. |4^, one can immediately 
see that both rephasing and non-rephasing parts of the 
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Figure 6: The four Feynman diagrams graphically illustrating 
the Liouville space pathways for a two-level system within 
the rotating wave approximation, g and e denote ground and 
excited state, respectively. 



real part 2D spectra are of a pronounced elliptical shape. 
However, the major axis of the ellipse in case of the 
rephasing part is oriented along the diagonal, whereas 
it is oriented along the anti-diagonal in case of the non- 
rephasing part. Each part on its own exhibits oscillations 
of its amplitude without significant line-shape modula- 
tions. The amplitude oscillations, however, are out of 
phase by tt, i.e. the rephasing part has its maxima at t2- 
times of 200, 450 and 650 fs, whereas the non-rephasing 
part has its maximal amplitude at t2 = 100, 300, 550, 
and 800 fs. Therefore, for ta-times of 200, 450, and 650 fs 
the major contribution to the total signal stems from the 
rephasing part, whereas for ta-times of 100, 300, 550, and 
800 fs the non-rephasing part contributes more strongly 
to the total spectrum. In Fig.[7|we present a quantitative 
evaluation of the respective contributions. 

An analogous situation is encountered in the imagi- 
nary part as shown in Fig. [lb. The nodal line of the 
rephasing part is oriented along the diagonal, whereas 
the one of the non-rephasing part is oriented along the 
anti-diagonal. As a consequence of the out-of-phase os- 
cillations of the amplitudes of the two contributions, the 
nodal line of the total signal experiences a periodic mod- 
ulation of its orientation. Similar to the real part, the 
rephasing and non-rephasing imaginary peak shapes do 
not change significantly with increasing ta-time, it is only 
their relative contribution to the total signal that varies 
with ^2- From inspection of the dispersive part it be- 
comes most obvious how the phase-twisted line-shapes 
in the rephasing and non-rephasing signal parts cancel 
upon summation to yield a purely dispersive imaginary 
part. In the accompanying paper II we explain the oscil- 
lations of the rephasing, nonrephasing, and sum spectra 
based on an expansion of the frequency-frequency corre- 
lation function in terms of Huang-Rhys factors. 

The dissection into rephasing and non-rephasing parts 
was also performed for the amplitude and phase spectra 
(Fig.O. For these representations, however, simple addi- 
tion of the two contributions to yield the total spectrum 
is not possible. Nevertheless, the oscillations in the to- 
tal spectra can be assigned to oscillating features in the 
rephasing and non-rephasing signal parts. In the ampli- 



tude part 2D spectra (Fig. [5^), we observe a stronger 
contribution of the rephasing part for ta-delays of 200, 
450; and 650 fs, and a reverse situation (i.e. an inten- 
sity gain of the non-rephasing contribution) for ta-delays 
of 100, 300, 550, and 800 fs. Therefore, the total sig- 
nal is elliptical and line-narrowed along the anti-diagonal 
in the former case, and more circular in the latter. In 
the phase spectra (Fig. [SJd), the lines of constant phase 
are oriented along the diagonal in the rephasing part and 
along the anti-diagonal in the non-rephasing part. This 
behavior reflects the different sign of the phase accumu- 
lated during the two time-intervals ti and according 
to the different coherence evolution in the rephasing and 
non-rephasing Feynman diagrams (59| . Depending on 
the relative strength of the rephasing and non-rephasing 
contributions, lines of constant phase in the total phase 
spectrum are either oriented along the diagonal (i2=200, 
450, and 650 fs), or they rotate down to become ap- 
proximately parallel to the wi-axis (^2=100, 300, 550, 
and 800 fs). In the next section we discuss relations of 
the line-shape evolutions in 2D electronic spectra to the 
frequency-frequency correlation function. 



C. Frequency- frequency correlation function 

For electronic spectra in the condensed phase, where 
transition energy fluctuations are non-Markovian in na- 
ture, homogeneous and inhomogenous timescales are not 
readily separable. Therefore, such systems are usually 
described by a frequency-frequency correlation function 
(FFCF) M{t), defined as 



M{t) 



{SuJeg it)duJeg{0)) 



(2) 



M{t) is a normalized ensemble averaged product of 
transition frequency fiuctuations {5ujeg) separated by 
time t. Since the correlation function in principle carries 
all relevant information on intramolecular and system- 
bath dynamics, considerable attention has been drawn 
over the last decade on how to sainple the FFCF by ex- 
perimentally feasible methods 0, [fiO, HH- Once M(t) 
is known, all linear and non-linear spectroscopic signals 
can be calculated via the line-shape function g[t) with 
knowledge of the reorganization energy A and the tem- 
perature dependent coupling parameter, which reduces 
at high temperatures. We use such complete 



to 



2\kBT 
h 



information about the transition frequency fluctuations 
due to intramolecular and system bath coupling to de- 
rive a detailed line shape theory in Paper II [3a|- Here, 
we concentrate on the intramolecular part of the M{t)- 
function which is responsible for the oscillatory dynamics 
of the line shapes. 

Although signal detection of one-dimensional FWM 
signals is relatively straightforward, these methods face 
an inherent trade-off between time- and frequency- 
resolution. This obstacle is evaded in 2D-ES, where 
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model calculations and experiments in the IR and VIS 
spectral region revealed correlations between line-shape 
evolutions and the frequency-frequency correlation func- 
tion [5i,[5i,|6iS|- 

Two quantities extracted from the real part of 2D elec- 
tronic spectra were shown to be directly proportional to 
the FFCF: the ellipticity of the peak shape [66] and the 
center line slope [67, 68]. The ellip ticity of the peak as a 
function of t2 is defined as [6^ [6g 



Eit2) 



(3) 



where D and A denote the diagonal and anti-diagonal 
peak widths, respectively. In Fig. [7^ we plot this ratio 
for < t2 < 800 fs, with the diagonal and anti-diagonal 
peak widths evaluated at ^-height. Oscillations with a 
period of 240 fs are clearly recovered. In our previous 
work we have shown that these oscillations arise mainly 
from modulations of the anti-diagonal width, whereas the 
diagonal width stays nearly constant over the time-range 
explored in our experiments [20]. The center line slope 
(CLS) is extracted by taking cuts parallel to the wa-axis 
and finding the maximum of the peak. The resulting 
data points as a function of wi are fit to a first order 
polynomial. The slope of this line is the center line slope 
(Fig. [7^) . Note that in the work by Kwak et al. [63| cuts 
are taken parallel to the wi-axis and the CLS is defined as 
the inverse of the slope of the line connecting the maxima 
of the slices. Both the ellipticity and the CLS can vary 
between 1 and and directly reflect the t2-dependent 
part of the FFCF . The limiting value of 1 can be 

reached for i2 = and in the absence of a homogeneous 
component. The deviation from the initial value of 1 can 
thereby be related to the magnitude of the homogeneous 
component, as was shown in Ref. [67'|. Relaxation and 
spectral diffusion processes during t2 lead to a decrease 
of the ellipticity and the CLS, until in the long time limit 
both quantities approach the value of zero, indicative of 
equal widths along the diagonal and anti-diagonal direc- 
tions and a slope of the center line that is parallel to 
the wi-axis. In a previous work it was argued that the 
CLS has the advantage over the ellipticity of being inde- 
pendent of factors influencing the line-shape (e.g. flnite 
pulse duration, apodization, etc.) [67| . As can be seen in 
Fig. [7^ we do not observe significant differences between 
these two quantities. Closely related to the ellipticity of 
the peak shape is the eccentricity, which is deflned as 

and which has been proposed in earlier 
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IF 



works as a measure for the FFCF [69|. Even though the 
eccentricity takes the same limiting values as the elliptic- 
ity, i.e. 1 if D>>A and if D=A (which is the case for 
a circle), the two quantities differ for intermediate cases 
and therefore cannot be compared directly. 

As discussed above, a periodic modulation of the 
rephasing and non-rephasing contribution to the total 
2D spectrum induces the observed changes in the real 
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Figure 7: Oscillating features extracted from the 2D spectra 
of PERY in toluene, (a) Ellipticity (filled circles) and center 
line slope (open circles) of the peak in the real 2D spectra, (b) 
Relative amplitudes of the rephasing (filled circles) and the 
non-rephasing (open circles) part of the real 2D spectra, (c) 
Inhomogeneity index extracted from the real part 2D spectra 
as defined in Eq. U) (d) Slope of the nodal line separating the 
positive and negative feature in the imaginary 2D spectra, (e) 
Slope of the = phase line in the phase 2D spectra. 
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and imaginary parts. In Fig. [7)d we plot the relative am- 
plitudes of the rephasing and non-rephasing part of the 
real 2D spectra. We clearly observe out-of-phase oscilla- 
tions of these two contributions, with the rephasing part 
exhibiting its maxima at i2-times of 200, 450, and 650 fs 
and its minima at i2 = 100, 300, 550, and 800 fs, and the 
non-rephasing part showing the opposite behavior, i.e. 
having its maxima at i2-tiines corresponding to minima 
in the rephasing part and vice versa. Tokmakoff et al. 
[5^ [65| defined an inhomogeneity index /(t2) as 



An- A 



NR 



Af 



A 



NR 



(4) 



with Ajf being the relative amplitude of the rephasing, 
and A^fi the relative amplitude of the non-rephasing 
part. This quantity was shown to directly relate to M(t) 
via M{t2) — sin(^). From the inspection of the in- 
homogeneity index (Fig. ^}p) one can notice, that the 
non-rephasing part exceeds the rephasing part in signal 
strength around i2-delays of 300 and 550 fs. For these 
delays the inhomogeneity index adopts negative values, 
which is consistent with negative peak shift values ob- 
served in Fig. [5^ (cf. Appendix). 

In the imaginary part of 2D spectra it is the slope of 
the nodal line separating the positive and negative con- 
tribution that is related to M{t) [64|. We extract this 
slope from our 2D spectra by finding the zero crossing be- 
tween the positive and negative contribution in the range 
-18700 cm-i < < -19200 cm"! and fitting of the re- 
sulting data points with a first order polynomial. The 
slope of the nodal line oscillates with a period of 240 fs 
(Fig. [TJi) . Note that also the slope takes negative values 
around t2 = 300 and 550 fs. This agrees well with the ob- 
servation of negative values of the inhomogeneity index 
(Fig. [71:) and the peak shift (Fig. [5^) . In case of a three- 
level system it is the real part that features two peaks 
of opposite sign, with the positive one stemming from 
ground state bleaching and stimulated emission and the 
negative feature arising from excited state absorption. In 
this case the slope of the nodal line is extracted from the 
real part spectra [H, [H, [t^ ■ 

Extracting the frequency-frequency correlation func- 
tion from the phase spectra (Fig. [SJd) is done by evaluat- 
ing the slope of lines of constant phase [5l,[6l|. In Fig. [7^ 
we plot this slope for = evaluated in the frequency 
range — 18500cm^^ < wi < —19250 cm^^. Again, the os- 
cillating pattern exhibits a period of 240 fs with maxima 
at t2-delays of 200, 450, and 650 fs and minima around 
t2 =100, 300, 550, and 800 fs. Around 300 and 550 fs the 
slope of lines of constant phase adopts negative values 
indicative of a tilt towards the anti-diagonal axis. 



V. CONCLUSION 

Intramolecular and system-bath (solvation) dynamics 
can be studied either in time-domain by the frequency- 



frequency correlation function M(t), or, equivalently, in 
frequency-domain by the spectral density p{u}). In this 
contribution we have concentrated on time-domain mea- 
sures, by presenting different metrics that can be ex- 
tracted from two-dimensional electronic spectra and that 
follow the functional form of the frequency- frequency cor- 
relation function. These metrics are the ellipticity, the 
center line slope, and the inhomogeneity index extracted 
from the real part spectra, the slope of the nodal line 
in the imaginary part spectra, and the slope of lines of 
constant phase in the phase spectra. Peak shift curves of 
PERY measured in solvents of different polarity (nonpo- 
lar, polar, and weakly polar) featured similar timescales 
[3J|, thus pointing to weak solute-solvent interactions. 
The observed oscillations in the frequency-frequency cor- 
relation function can thus be assumed to be dominated 
by intramolecular dynamics. The observed beating pat- 
tern stems from the fact that during t2 the density ma- 
trix does not only evolve in a population state (ground 
or excited state), but vibrational coherences can be ex- 
cited as well. These coherences evolve with a frequency 
that is given by the difference in frequency of the two 
states involved. In theoretical works on vibrational ef- 
fects in 2D-ES, the vibrational sub-levels could be re- 
solved as a consequence of the customized input param- 
eters [13, [3^ • Under experimentally feasible conditions 
effects as finite pulse durations and finite temperature 
erode the multipeak structure. In such a situation, vi- 
brational dynamics can not be followed by the evolu- 
tion of diagonal and cross peaks, but by the evolution 
of line-shapes and signal intensities, as has been done 
in the present contribution. In our previous work (20| 
we performed a fit to the experimental two-dimensional 
spectra by including two overdamped solvent modes and 
two underdamped intramolecular modes. In the second 
paper of this series we investigate the influence of vibra- 
tional coupling on 2D electronic spectra more thoroughly 
by employing nonlinear response function theory and ex- 
panding the frequency-frequency correlation function in 
terms of Huang-Rhys factors. Furthermore, we extend 
the theory also to the case of high frequency modes in 
which case individual diagonal and cross peaks can be 
distinguished. 



Acknowledgements. This work was supported 
by the Austrian Science Foundation (FWF) within 
the projects No. P18233 and F016/18 Advanced 
Light Sources (ADLIS). A.N. and J.S. thank the Aus- 
trian Academy of Sciences for partial financial support 
by the Doctoral Scholarship Programs (DOCfFORTE 
and DOC). T. M. acknowledges the kind support by 
the Czech Science Foundation through grant GACR 
202/07/P278 and by the Ministry of Education, Youth, 
and Sports of the Czech Republic through research plan 
MSM0021620835. J. H. gratefully acknowledges support 
by the Lise Meitner project M1080-N16. The quantum- 
chemical calculations were performed in part on the 



11 



(a) (b) (c) 




t2 / fs t2 / fs t2 / fs 

Figure 8: One-dimensional four-wave mixing signals of PERY in toluene. Top and bottom panels share the same abscissa, color 
scales are shown as insets in the bottom panels, (a) Three-pulse photon echo peak-shift trace (top) and the corresponding 
frequency integrated three- pulse photon echo signal as a function of delay ti and t2. (b) Frequency integrated (top) and 
frequency resolved (bottom) transient grating signal, (c) Frequency integrated (top) and frequency resolved (bottom) pump- 
probe signal. 



Schrodinger III cluster at the University of Vienna. 



Supplemental Material. The Cartesian 

B3LYP/SV(P) coordinates, total electronic energy, 
harmonic vibrational spectrum, Huang-Rhys factors, 
and a visualisation of the most dominant vibrational 
modes of PERY can be found in the supplemental 
material. 



Appendix A: Comparison to other four-wave mixing 
signals 

Heterodyned two-dimensional electronic spectroscopy 
gathers the maximum amount of information that can 
be inferred by any third-order nonlinear technique. Nev- 
ertheless it is instructive to cross-check our findings 
against results that are obtained with alternative detec- 
tion schemes, shown in Fig. [5] From the experimental 
point of view, these signals may be differentiated by the 
required measuring (scanning) time and set-up complex- 
ity, while, on the other hand, all of them are related to 
the functional form of M{t) UM,^- 

The three-pulse photon echo peak-shift (3PEPS) 
method is a variant of homodyned photon echo spec- 
troscopy, in which the time-integrated nonlinear signal 



S3PE{tl,t2)^ / \P'^^\h,t2,h)fdts (Al) 
JQ 



is recorded in a wavevector architecture equivalent to 2D- 
ES. In each step of a 3PEPS sequence, delay ti is scanned 
at a fixed delay t2, and the position of the (time- and 
frequency-integrated) signal maximum evaluated along 
the ii-axis. The magnitude of this peak-shift (devia- 
tion from <i 0) is recorded for a range of f2-delays, 
to give the peak-shift decay. In essence, a finite pho- 
ton echo peak-shift value indicates the system's ability to 
rephase a macroscopic coherence (i.e. to cause an echo) 
after having spent a certain time (^2) in an intermediate 
state. Fig. [5^ shows the photon echo peak-shift decay 
for PERY in toluene, along with a two-dimensional plot 
of S3PE{ti,t2)- Both data representations clearly reveal 
how the photon echo signal maximum oscillates (along 
ti) as a function of t2, the oscillation being observable 
for i2-delays well beyond one picosecond. In the context 
of our analysis of the 2D electronic spectra of PERY, it is 
interesting to note that the peak-shift acquires negative 
values around the minima of the peak-shift trace (located 
at i2-delays around 100, 350, and 600 fs). In other words, 
the integrated photon echo (recorded as a function of ti) 
may peak at negative ii-delays, if non-rephasing signal 
contributions dominate over rephasing contributions to 
the total signal. We point out that this finding is consis- 
teirt with the negative value for the inhomogeneity index 
and a negative slope of the nodal line in the dispersive 
2D signal part, observed at the corresponding i2-delays. 

For the sake of completeness. Fig. [8}d and c show the 
transient grating (TG) and the pump-probe (PP) signals 
of PERY, in both frequency integrated and frequency re- 
solved representation. The former method records S^pe 
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as a function of t2 for delay ti set to zero (descriptively, 
the first two excitation pulses form a spatial transmis- 
sion grating from which the third pulse is scattered off). 
Note that unlike 3PEPS, TG measurements are sensitive 
to population decay during i2- Time delays are scanned 
in a similar way in PP measurements (i.e. scanning t2 
for ti — 0), however, only two pulses are involved in the 
experiment, the signal being recorded along the direction 
of the second pulse. Since the PP signal is generated by 
two interactions with the first pulse (pump) , followed by 
one interaction with the second pulse (probe), ti is intrin- 
sically zero. Overall, also our TG and PP data-sets do 



well align with previously published data [1, [s^l and our 
2D electronic spectra. In both types of signals, strong 
intensity oscillations with a period of 240 fs (140 cm~^) 
can be discerned. The slow decay of the TG signal con- 
firms the absence of any fast population decay channels. 
In line with the discussion above, the PP signal is domi- 
nated by (positive) stimulated emission and ground state 
bleaching contributions, indicating excited state absorp- 
tion effects to be negligible on the picosecond timescale 
investigated here. 
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